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We construct an explicit model to describe fermions confined on a four dimensional brane em- 
bedded in a five dimensional anti-de Sitter spacetime. We extend previous works to accommodate 
massive bound states on the brane and exhibit the transverse structure of the fermionic fields. We 
estimate analytically and calculate numerically the fermion mass spectrum on the brane, which we 
show to be discrete. The confinement life-time of the bound states is evaluated, and it is shown that 
existing constraints can be made compatible with the existence of massive fermions trapped on the 
brane for durations much longer than the age of the Universe. 
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I. INTRODUCTION 



The idea that our universe may be a hypermembrane in a five dimensional spacetime has received some attention in 
the last few years after it was realized that gravity could be localized on a three-brane embedded in an anti-de Sitter 
spacetime jl] . Since then, much work has been done in a cosmological context Q and there is hope that a consistent 
(i.e., mathematically self-contained and observationally satisfying) high dimensional model might soon be formulated. 
For instance, it has been proposed that a model || based on such ideas could present itself as an alternative to the 
inflationary paradigm, although for the time being the controversy as to whether or not such a model might have 
\ anything in common with our Universe is still going on Q . 

The idea is not new however, but has evolved from the standard Kaluza-Klein approach to that of particle local- 
O i| ization on a higher dimensional defect (|||. In particular, it has been shown that massless bulk scalars and gravitons 
share the property to have a zero mode localized on the brane in the Randall-Sundrum model. Various mecha- 
'-^ ' nisms [|| have been invoked according to which it would be possible to confine massless gauge bosons on a brane, so 
. that there is hope to achieve a reasonable model including all the known interactions in a purely four dimensional 
J/Tj ' effective model. 

A mechanism permitting localization of massless fermions on a domain wall was described in Refs. I^Jj. However, 
■ although appealing this mechanism might be, it should be emphasized that actual fermions, as seen on an everyday 
basis in whatever particle physics experiment, are massive, so that a realistic fermionic matter model on the brane must 
accommodate for such a mass. The question of localization of massive fermions on the brane thus arises naturally, and 
it is the purpose of this work to provide the transverse brane and fermionic structure that leads to this localization. 
Up to now, fermions have been confined under the restricting hypothesis that the brane self gravity was negligible [jul 



or that it was embedded in a Minkowski spacetime with one 11 1 or two [|12| transverse dimensions (see also |13[ for 
the localization of fermion on a string-like defect in five dimension). 

Our goal is to transpose the original work of Ref. Jl2| to the brane context. For that purpose, we realize the brane 
as a domain wall. Such domain wall configurations in anti-de Sitter space have already been studied jl4|,|l5| . We will 
assume that five dimensional fermions are Yukawa- coupled to the domain wall forming Higgs field, as in the usual 
case of cosmic strings. In this respect, our work somehow extends Ref. |lC[ ], where the mass term was put by hand, 
and Ref. Jll| where the gravity of the wall was neglected. 

We start, in the following section, by recalling the domain wall confi gura tion of a Higgs field in a five dimensional 



anti-de Sitter spacetime and discuss briefly its properties. In section III, we describe the dynamical equations of 
fermions coupled to this domain wall in order to show that they obey a Schrodinger-like equation with an effective 
potential which can trap massive modes on the wall. The asymptotic structure, i.e., deep in the bulk (far from the 
brane), is not Minkowski space, so that the effective potential felt by the fermions possesses a local minimum at the 
brane location, but no global minimum, as first pointed out in Ref. E(j ]. As a consequence, the bound states are 
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metastable and fermions can tunnel to the bulk. 

We then provide an analytical approximation of the effective potential, thanks to which we compute analytically, 
in section IV, the mass spectrum of the fermions trapped on the brane. We obtain the mass of the heaviest fermion 
that can live on the brane and estimate its tunneling rate. This result is compared to a full numerical integration, 
performed in section |v|. In a last section, we investigate the parameter space and, after having compared our results 
to previous ones, we conclude that there exists a wide region in the parameter space for which the fermion masses 
can be made arbitrary low, i.e., comparable to the observed small values (with respect to the brane characteristic 
energy scale), while their confinement life-time can be made much larger than the age of the Universe. Such models 
can therefore be made viable as describing realistic matter on the brane. 

II. MEMBRANE CONFIGURATION IN ADS 5 

We consider the action for a real scalar field $ coupled to gravity in a five dimensional spacetime 

S = J ^{R-2K) + ^g AB d A <S>d B 'S>-V{<$>) V5d 5 x = y"v?d 5 x[£ grav +jC a + U] (1) 

where qab is the five dimensional metric with signature (+,—,—,—,—), R its Ricci scalar, A the five dimensional 
cosmological constant and k 2 = 67r 2 G 5 , G 5 being the five dimensional gravity constant. Capital Latin indices A,B... 
run from to 4. The potential of the scalar field $ is chosen to allow for topological membrane (domain wall like) 
configurations, 

U($) = ^($ 2 - f7 2 ) 2 , (2) 

where A is a coupling constant and r\ — (|<I>|) is the magnitude of the scalar field vacuum expectation values (VEV)Q. 
Motivated by the brane picture, we choose the metric of the bulk spacetime to be of the warped static form 

ds 2 = g AB dx A dx B = -dy 2 + e^^^da^d^ = -dy 2 + g^dx^dx", (3) 

where r]^ is the four dimensional Minkowski metric of signature (+, — , — , — ), and y the coordinate along the extra- 
dimension. Greek indices /i, v . . . run from to 3. 

With this metric ansatz, the Einstein tensor components reduce to 

G^ = - gflu (6a' 2 -3a"), G yy = -6a' 2 , (4) 

where a prime denotes differentiation with respect to y. The non- vanishing components of the matter stress-energy 
tensor 

Tab = 2 ^ A % ~ 9abC^, (5) 



are given by 



T M u = \g^' 2 + 2V), Tyy = i($' 2 - 2V). «ii 



It follows that the five dimensional Einstein equations 

Gab + kg AB = k\Tab (7) 

can be cast in the form 



Note, that, because of the unusual number of spacetime dimensions, the fields have dimensions given by [7?] = M 2 , [$] =M 3/2 , 
[A] = M 2 , [A] = M"\ [77] = M 3/2 and [k b ] = M" 3/2 (M being a unit of mass). 
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= $' 2 (8) 



K 2 



6(T ' 2 = *S. ($' 2 - 2V) - A, (9) 

while the Klein-Gordon equation takes the form 

$"-4<7'$' = ^. (10) 

Eqs. (§-f[o|) is a set of three differential equations for two independent variables ($ and a). Indeed, as can easily be 
checked, the Klein-Gordon equation stems from the Einstein equations provided $' ^ 0. To study the domain wall 
configuration, we choose to solve the first Einstein equation (||) together with the Klein-Gordon equation ((Tof). 

This system of equations must be supplemented with boundary conditions. By definition of the topological defect 
like configuration, we require that the Higgs field vanishes on the membrane itself, i.e., $ = for y = 0, while 
it recovers its VEV in the bulk, so that lim^-too <£> = ±77. Note that the sign choice made here is arbitrary and 
corresponds to the so-called kink solution ; the opposite choice (i.e. \hn y ^± 00 $ = ^pq) would lead to an anti-kink 
whose physical properties, as far as we are concerned, are exactly equivalent. As for the metric function a, it stems 
from the requirement that one wants to recover anti-de Sitter asymptotically, so that one demands that a' tends to 
a constant for y — > ±00. This constant can be determined using Eq. (g), so that lim^-too a' — ±y- A/6. Note that 
as y changes sign at the brane location, there is no choice for the sign of the function a in this case. Note also that, 
as is well known, the static hypothesis implies that the bulk cosmological constant A must be negative, and therefore 
the five dimensional spacetime to be anti-de Sitter. 

With the convenient dimensionless rescaled variables 

g = yV\A\, H = -, S=^, (11) 
77 dg 



the dynamical equations read 



3# 2 , (12) 



H - ASH = Aj3H{H 2 - 1), (13) 

where a dot refers to a derivative with respect to g and the two dimensionless (positive) parameters a and (3 are 
defined by 

(") 

These parameters are not independent since, for an arbitrary value of (3 say, there is only one value of a for which 
the boundary condition 5(0) = 0, or equivalently lim e _ + _ 00 S(g) = — is satisfied. This stems from the fact that 

Eq. (^2|) is a first order equation in 5, so that only one boundary condition is freely adjustable, and we choose it 
to be at g — * +00. Once this choice is made, the value of 5 on the brane is completely determined, and unless the 
parameters are given the correct values, it does not vanish. As the solution must be symmetric with respect to the 
extra dimension coordinate y, one must tune the parameters in order to have a meaningfull solution (i.e. for which the 
metric and its first derivative are continuous at g = 0). This is reminiscent of the relation that should hold between 
the brane and bulk cosmological constants pl| . 
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FIG. 1. The rescaled Higgs field amplitude H (full line for /3 = 0.01 and dashed line for (3 = 0.1) and warp factor derivative 
S (dotted line for f3 = 0.01 and dot-dashed line for f3 — 0.1) as functions of the rescaled extra dimension coordinate g. 

Eqs. (|l2| - |l3| ) have been solved numerically with the relevant boundary conditions. The field profiles are depicted 
on Fig. |I]for two arbitrary values of the parameter (3. 




P 

FIG. 2. Relation between the dimensionless constants a and /3. The thick curve represent the result of the numerical 
integration, while the (hardly distinguishable) thin curve is the best analytical fit. 

The relation between the parameters a and (3 such that the metric is regular at the brane location is depicted on 
Fig. 0. As can be seen on the figure, it consists essentially in two power laws. For small values of /?, one finds roughly 
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a ~ 1/(3, which becomes exact in the limit /3 — ► 0, while for large values of /3, one gets a ~ |/3 1 / 2 , which is, again, 
exact in the limit /? — >• oo. We were able to find the best fit 



(15) 



which, as can be seen on the figure, is almost exact everywhere. This translates into the relation 



|A| = qAt? 




(16) 



between the 5 dimensional cosmological constant and the microscopic parameters, which corresponds to the usual 
relations between brane and bulk cosmological constants. 

III. LOCALIZATION OF FERMIONS ON THE WALL 

This section is devoted to the description of the Dirac equation in five dimensions with the domain wall configuration 
obtained in the previous section. 

The minimal representation of spinors in five dimensions can be chosen to be four dimensional Jl6| . The five dimen- 
sional Clifford algebra can then be constructed from the usual four dimensional one by adding the 75 matrix to close 
the algebra. If 7^ design the usual four dimensional Dirac matrices in Minkowski space in the chiral representation, 
the Dirac matrices in five dimensional Minkowski space, T A , are 

r» = y, r 4 = -i 75 (17) 

and they satisfied the usual Clifford algebra 

{T A ,T B } = 2 V AB (18) 
where T] stands for the five dimensional Minkowski metric. Since in this representation the Dirac matrices satisfy 

r 5 = ir°...r 4 = id (19) 

the five dimensional spinors have neither Weyl nor Majorana represention. It follows that the Dirac Lagrangian in 
five dimension for fermions coupled to the Higgs domain wall is necessary of the form 



= {WT A V A ^ - 5f *$*) = y/gV [ie a{v) T^V^ + iT 4 d v - g F $ 
where the Lorentz covariant derivative with spin connection is fr?j| 



(20) 



(21) 



We emphasize that the sign of the coupling g F of the spinor ^ to the Higgs field is arbitrary and represents a coupling 
either to kink or to anti-kink domain wall. For definiteness, we shall consider in what follows only the case of a kink 
coupling, and thus assume without lack of generality that g F > 00. 

The variation of the Lagrangian (^0|) leads to the equation of motion of the spinor field, namely the Dirac equation 
in five dimensional anti-dc Sitter space for a fermionic field coupled to a Higgs field, 



iT 4 (d y - 2a' (y)) + ze^r^ - 5f $ 



* = 0. 



(22) 



This equation involves the matrix 75 (through 1^4) and it is thus convenient to split the four dimensional right- and 
left-handed components of the five dimensional spinor and to separate the variables as 



"Note also that the dimensions are given by [*] = M 2 and [g F ] = M 



-1/2 



5 



*(x»,y) 



1 + 75 



Uk{v) 



1-75 



(23) 



where ip(x^) is a four dimensional Dirac spinor, while U R {y) and U\Xv) are yet-underdetermincd functions of y. 
In what follows, we want the five dimensional Dirac equation to yield an effective four dimensional massive Dirac 
equation, with an effective mass m (energy eigenvalue of the bound state). Such a requirement implies that 



ij^d^ip — mip, 

or, equivalently, in terms of the right- and left- handed components 

i^d^ipn = mip h , i^d^ipt = mtpn, 
where the right- and left-handed components of the four dimensional spinor are defined as 



1 



75 



and 



1-75 



(24) 
(25) 

(26) 



Contrary to the case studied in Ref. |Tc| ] , the mass m is not an arbitrary parameter and will be determined later. 

Choosing ?/>r and as the independent variables instead of \& and and inserting equation ( p5|) into the equation 
of motion ( p2| ) while using the splitting ansatz (^3|) yields the differential system for the two functions Wr/l(2/), 



[d v - 2a' (y) - g F <P]U R (y) = -me'^fe), 
[d y - 2a'(y)+g F ^>]U L (y) = me^U R (y). 



(27) 
(28) 



To simplify the notations, it is convenient to introduce the dimensionless rescaled bulk components of the fermions 



in terms of which the system fl2 



U{q) = 
takes the form 



|A|V4' 



9 S - 



d„ 



Ml 



with the dimensionless rescaled mass and coupling constant 

m 



and 7 F 



9fV 



(29) 

(30) 
(31) 

(32) 



Let us first concentrate on the special case \i = 0. The system (|30|-|3lD then consists in two decoupled differential 
equations and the zero mode states H] are recovered. Asymptotically, these functions behave as 



W R (p^±oo)^e( 7F+ ^) lel , 



U L (g -> ±oo) ~ e v 7p ^ 



I lei 



Thus, only the left-handed solution Ul may remain bounded @||, and yet provided 

f 



7 F > 



2V6' 



(33) 
(34) 

(35) 



Indeed, the right-handed zero modes could have been obtained by considering the coupling of fermions to the anti-kink 
Higgs profile^. We thus recover the well-known fact that massless fermions must be single-handed in a brane model, 
contrary to the ordinary four dimensional field theory in which they simply can. 



3 The coupling with and anti-kink for which g F < would have yield right-handed solution with the constraint 7 F < — 



G 



Let us now focus on the more interesting massive case for which /i ^ 0. Then the system (|30jj3l|) can be decoupled 
by eliminating say. For that purpose, we differentiate ( j3l"l ) with respect to g and express OqUr using equation ( |30| ) 
and U-r. using equation ([HJ) again to get 



d 2 e -2S0 e 



Z4=0, 



da 



(36) 
(37) 



This system is strictly equivalent to the initial system (^-^Tj) since differentiating Eq. ( p7| ) and then using Eq. ( |3"6| ) 
gives back Eq. (|3(]). It is thus important to keep both equations. Note that the integration of the first equation 
( |36| ) will require two initial conditions but that Wr will then be completely determined and thus requires no extra 



constant of integration. As a consequence, it is sufficient to solve the second order equation (36) for Ul in order to 
fully determine the left- and right-handed bulk fermion profiles. 

Eq. (Bq) can be recast into a Schrodinger-like second order differential equation 



d 2 U h +to 2 (g)U h = 0, 



where the function to is defined by 



with the new function 



j 2 (g) = M 2 e 2 ^> + d e ( % H + ±s) (j F H + ±S 



U{g) = e- a ^U(g). 



(38) 



(39) 



(40) 



Our aim will now be to find the zero modes of this new equation; as previously discussed, they will be equivalent to 
the massive bound states we are looking for on the brane. 

In order for the fermions to be confined on the brane, the minimum of lo 2 needs to be negative to imply an 
exponential decrease of £/l in the bulk. This is essentially equivalent to the condition ( |35| ) that was obtained for the 
case of zero modes. We shall assume henceforth that this condition also holds for massive modes, i.e., that the value 
of 7 F necessary to bind massive fermions on the brane is at least that to bind massless ones. Indeed, Eq. ([39j) shows 
that the minimum of lo 2 can only be negative for large values of the parameter 7 F . However, since the first term 



of (39) increases exponentially at large distance from the brane, then, if [x ^ 0, lo will necessarily become positive. 



This will yield asymptotic radiative behaviors of the spinor bulk components. Physically, it can be interpreted as a 
tunneling of the fermions from the brane to the bulk 113] . On the other hand, on the brane, the Higgs field H and 
the derivative of the warp factor S vanish, so that u> 2 (0)is positive. As a result, the fermions can freely propagate in 
a tiny region around the brane, but certainly only for particular values of io 2 (and thus of \i) satisfying the boundary 
conditions with the surrounding exponential decreasing regions. The effective potential V a{{ — —lo 2 , depicted in Fig. |^, 
exhibits a local minimum on the brane and minima at infinity. The modes trapped on the brane are thus expected 
to have discrete masses /i on the brane and non-zero probability of tunneling into the bulk. 
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FIG. 3. The effective potential V ctf = — uj 2 felt by Wl- We chose the parameters 7 F = 100 and fi = 0.5j F (an illustrative 
value, not necessarily leading to a bound state) in the Higgs and gravity background similar to those of Fig. hi and we have 
assumed, for definiteness, the values a = 1.63, j3 = 1.23. The fermions are trapped on the brane where V off is negative and have 
a non-zero probability of tunneling into the bulk due to combined effects of Higgs and gravity which produce a finite potential 
barrier. 



IV. ANALYTIC ESTIMATE OF THE MASS SPECTRUM AND OF THE TUNNELING RATE 



Since the Higgs and warp factor profiles are not known analytically, it is a priori impossible to solve Eq. (|38 
analytically. Nevertheless, since the fermions are expected to be trapped in a neighborhood of the brane, we can look 
for series solutions in g. 

In § IV A, we give an approximat ion of the effective potential V ct{ which will then be us ed to determine the bound 

We end this section by determining in § IV C the tunneling rate in this 
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states and the mass spectrum in 
approximation. The validity of this approach is difficult to assess and it will be justified a posteriori on the ground of 
a full numerical integration of the system in the following section. In the whole section, it is assumed that fi > 0. 



A. Approximation of the effective potential 



In a neighborhood of the brane we can expand the Higgs and warp factor profiles as 

H(g) = H ig + H 3 g 3 +0(g 4 ), 
S(g)^S 1 g + S 3 g 3 + 0(g i ), 



(41) 
(42) 



where both the constant and quadratic terms vanish for symmetry reasons. To simplify the analysis, we shall make 
use of the equations of motion in the form 



S=-F(S,H), 

and 

H = ±J-F(S,H), 
V a 

where the function F(S, H) is defined as 

F(S, H) = 6S 2 + a(3{H 2 - if - 1, 



(43) 



(44) 



(45) 



(recall that A < 0). 

Plugging the expansions ( pd] ) and (^) into Eqs. ( f43| ) and (H), it follows that the function F(S, H) can be expanded 
up to third order as 



F(S, E)~ a p-1+ (6Sf - 2aPHl) g 2 + 0{g 4 ). 
Inserting this expression back into the equations of motion yields the three coefficients 



Si 



-H 



i ) 



H 3 - -H, 



;Hl (3 



in terms of the coefficient Hi 



Hi = d e H\ e ^ = \j-(ap-l). 



(46) 

(47) 
(48) 
(49) 

(50) 
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Then, the frequency u> 2 can be expanded as an harmonic oscillator potential 



^{ Q ) = ul{Q) + o{ e i ) 



ith 



(51) 



where u!q and fl are given by 



2 _ 2 
CJ = /i 



(52) 
(53) 



The function w 2 is well approximated by only near the brane and the expansion ( pl| ) is no longer valid at large 
distance where the exponential term dominates [see Eq. (39)]. Once the fixed asymptotic values of the Higgs and 
warp factor are reached, the frequency (J39l) behaves as 



with 



7 F + 



2%/6 



(54) 



The analytical estimate of the function uj 2 is thus obtained by matching the two limiting asymptotic behaviors (|51 
and (M), respectively lo\ closes to the brane and uj 2 ^ far from it, as 



co 2 (p) 



< \P\ < Pu 
<o \p\ > P a 



The dimensionless matching distance p m has to be solution of 



(55) 



(56) 



in order to get a continuous function. Note that, because of the symmetry on both sides of the wall, we can assume 
g m > without lack of generality. For large values of 7 F , and using Eq. d53) , we get 



leading to 



1 

Hi' 



(57) 



(58) 



As it turns out, the faster the asymptotic solution is reached, the better the approximation works. This is the case in 
particular for H\ > 1, 

The exact (numerically integrated) effective potential and its approximation are compared in Fig. ^. The global 
shapes are effectively the same, and in spite of uncertainties at intermediate regions due to this crude approximation, 
it is reasonable to expect the same fermion physical behaviors in both potentials. Note that for (cosmologically 
favored) higher value of a/3, the Higgs field and the warp factor reach more rapidly their asymptotic values leading 
thus to a better agreement between the two potentials, as can be seen on Fig. ||. 



B. Determination of the bound states 



Given the approximate frequency (k33) , the equation of motion for the left-handed bulk spinor component ti^ reduces 



to 



Q > Qn 



Q < Qn 



a 2 + /x 2 e 2| e |/V6_ (7F + 



[d 2 + lu 2 - Qg 2 ] U h - 0, 



2^/6 



Wl = 0, 



(59) 
(60) 



with the requirement that 11^ and its derivative are continuous at g m . Note also that we consider only the case g > 0, 
physics on both sides of the brane being completely symmetric under the transformation g — * —g. Let us consider 
the solutions in each region separately. 
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FIG. 4. The effective potentiai V off = —uj 2 and the approximate analytic effective potential (dashed curve) obtained from 
the matching of its two asymptotic expansions near the brane and at infinity. On the left panel, the parameters are the same as 
on Fig. ^, i.e., 7 F = 100 and ^ = 0.57 F . As the asymptotic solution is not reached very close to the brane, the approximation 
is rather poor and the two solutions are not in good agreement at intermediate regions. On the right panel, we use the new 
parameter values 7 F = 100 and /i = 0.67 F , obtained for a = 0.1865 and f3 = 53.62, a much better approximation is obtained 
around the potential barrier provided the Higgs field and warp factor reach their vacuum value rapidely, i.e., for larger value of 
a(3. As on Fig. [], the parameters are chosen to illustrate the point and do not necessarily correspond to existing bound states. 



• 8 > Q m '■ we introduce the new variable 

z = V6fie lsUV ^, 

in terms of which Eq. ( |59| ) reduces to a standard Bessel differential equation 

W L = 0, 



dz 2 z dz 



the order of which being given by 



V&1 F + \- 



(61) 



(62) 



(63) 



Since uo 2 is positive at infinity, the asymptotic form of the solution is necessary radiative, as was already pointed 
out in the previous section. The most general solution of the Bessel equation (62) is a linear superposition of 



Hankel functions. Since we are interested only in ingoing waves in order to study a tunneling process, the most 
general solution takes the form 

Z40) = BH<f\z), (64) 
where Hf~\z) is the Hankel function of the first kind, propagating towards the brane at infinity fl9[| 

(65) 



^ (1) (z^oo) 



2 e i(z-£7r/2-7r/4) 



and B is an arbitrary complex constant. 
• Q < Q m '■ performing operations similar to those of the previous case, we cast Eq. ( |60|) on the form 

" d 2 (1 



Z4 = 0, 

in which we have introduced the new variable and parameter 



(66) 
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1 J 1 

s = (4n)' e , a ee 2_ (67) 

The general solutions of Eq. ( p6| ) are the parabolic cylinder functions, namely U(a,x) and y(a, x), of which 
lhj{x) can be expressed as linear superpositions. In the limit x |a|, these solutions scale as p5[ 

l7(o,a?) ~e-* x V°-', and V(o,x) ~ J^e^x?-*. (68) 

Since we are interested in confined fermion states on the brane, only the exponentially decreasing function is 
relevant, so that the general solution near the brane reads 

Ul{x) = AU(a,x), (69) 

where A is a complex integration constant. 

The general solution Ul(q) for all g is obtained by matching the two different solutions at g = g m . Since g m 
corresponds to the maximum positive value of the effective potential [see Fig. [I| , it is reasonable to consider that the 
Hankel function at that point can be expanded around small values of their argument with respect to their order [ fl9| , 
i.e., 

H { l\z m )~--T{l)z- i , (70) 

while the parabolic cylinder functions can be taken in their large argument asymptotic limit (|68| ) . This is the same kind 
of approximation as that made to derive the effective potential. Physically the initial conditions on the brane, i.e., 
Wl(0) and d g U^\o, are chosen in such a way that the asymptotic exponentially growing function V(a,x) contribution 
is everywhere negligible. Once these initial conditions are fixed, they fully determine the solution on the other side of 
the brane, i.e., for x < 0. The asymptotic expansion (pq) can be analytically extended to — |x| = \x\e™ and yields p]J| 

U(a, -\x\) ~ e-iW 2 \x\- a -h- t7r ( a+ ^. (71) 

Thus, once Wl(0) and d e Uh\o are fixed, the matchings between Hf l \z m ) and U(a, — \x m \) on one side, and U(a, \x m \) 
on the other side fully determine the bulk component for all g. 

The last constraint comes from Eq. ( |37| ) determining the right-handed spinor bulk function. It is well defined if 
and only if both IA\^ and d g Uh are not singular. In fact, the derivative of the parabolic cylinder function U(a, x) is 
generally discontinuous at x = 0. With the help of the Wronskian of U(a,x) and U(a, — x) J2(]] 

dU(a,-x) dU(a,x) 2ir 

t/(fl ' X) ^x~ - u{a > - X) ^x~ = T{rT^y (72) 

we can construct the derivative discontinuity at x = 0. This is 

d£/(a,0-) dt7(a,0+) _ . + 3 r(| + §) 

dx dx r(i+a)' 1 ' 

where we have used the particular value p0| 

U(a, 0) = - ; ^ -. (74) 

2 f+lr(| + f) 

Imposing that the derivative of is continuous at x — results in imposing that the jump (|73) vanishes. This is the 
case if and only if a is solution of 



r(f + f) 



r(i + a) 

Since T is singular for negative integer arguments, this condition is satisfied only for 



0. (75) 
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where n is a positive integer. Note that —a — 1/2 cannot be odd since then the numerator of the Wronskian ( |75| ) will 
also be singular resulting in a finite derivative jump at x = 0. The condition (|7^ ) shows that the trapped fcrmions 
on the brane have necessarily discrete masses fi n which read, using the values of the parameters ([52]), (^3|) and (|67j), 



(4rc + 1)4/1 + 



6/3 24a(3 + a 2 Hf [(4n + l) 2 - 5] 



3 7f #i 



36 7f 2 



1 +-^[( 4 "+ 1 ) 2 + 1 ] 



(77) 



This mass spectrum is valid for n > since our derivation assumed that [i > 0. In the limiting case where 7 F 3> 1, it 
reduces to the much simpler form for the lowest masses 



(78) 



On the other hand, /i 2 cannot reach very large values since it is necessary to have a potential barrier in order to have 
bound states. From the expression of the effective potential (EJ), the barrier is found to disappear when 



^(^■fmax) ~ 0. 



(79) 

Again in the limit where 7 F 3> 1 , using the value ( j58| ) of £> m , one gets the maximum accessible reduced mass /i m ax for 
fi as 



rmax if c 

The maximum number of distinct massive states trapped on the brane can thus be estimated to be 



(80) 



Int 



1 



7f 



4iJi 



For the parameters chosen in Fig. |^, one obtains /imax ~ 0.687 F and there are n max = 11 massive modes trapped on 
the brane. 



C. Fermion tunneling rate 



Since the effective potential becomes negative at infinity, the massive modes trapped on the brane are subject only 
to a finite potential barrier. They are in a metastable state and can tunnel from the brane to the bulk. In this section 
we use our previous analytic solution to estimate the tunneling rate. Would this rate be too high, one would observe 
an effective violation of energy-momentum conservation on the brane, i.e., in four dimension, thereby contradicting 
observation. 

Our starting point is the analytic solution for the left-handed bulk function that was derived in the previous section 



e > Q n 



U l (q) = AU (a, V20 1A 
U L (g)=BH i i 1) 



(82) 
(83) 



The transmission factor can easily be derived from the matching conditions of the left-handed bulk fermion component 
at g = g m . First of all, Ul has to be continuous. Using the expansions ([38]) and ( |70| ) and the value ( |76| ) of the parameter 
a, permits to find the relation 



between the coefficients A and B. Making use of the expression (|58|) for g m yields 




(84) 



B in 



A T(t) 




exp 



12Vfi 



(85) 
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Assuming, as above, that 7 F ^> 1 (so that £ 3> 1) we can expand T(£) as 



(86) 



so that, using expressions (57), ( J5q ) and ( p3[ ) respectively of fi, p m , and we get that B/A is approximated by 

2 



B . fir 
A ~ V 2 V 



gl/4 n+1/2 

' F 



exp 



_V6 7f (lnH^_i+ L_ 
V M 2V6i/i 



(87) 



The transmission coefficient from the brane to the bulk associated with can thus be defined by 

Z4(£) 



T 



Z4(o) 



with Ul(£) — BH^\l), is evaluated at the turning point g = £ where the spinor bulk component begins to propagate 
freely. Using the behaviour (|74|) of the function at the origin, the ratio (87) and the properties Jl9] of the Hankel 
function 



H^(£) 

the transmission coefficient (|8S| ) reduces to 
f 



2/3 



r/3 



r(2/3) 



2V3r(l/2-n) ; „ +1/ 
3V«r(2/3) 1 7f p 



-V&Tf ( In 



2Md 



/' 



2V6i?i 



(89) 



(90) 



It follows, using the definition (40), that the probability for a trapped particle on the brane to tunnel to the bulk is 
given by 

(91) 



v = \r\ 2 = c 2e ^ e \t\ 2 . 

The characteristic time for a fermionic mode trapped on the brane can be roughly estimated by 

r b ~y b = g » 



(92) 

where y b represents the typical length, in the fifth dimension, felt by a particle on the brane. As can be seen on Fig. |^, 
the spinor bulk components are exponentially damped as soon as the effective potential becomes positive. Thus g b 
can be estimated by the solution of u> 2 (g b ) = so that, keeping in mind that \i < j F , 



The life-time r„ of a fermionic bound state on the brane labelled by n 



V 



can be estimated by 



3 7 / 6 r 2 (2/3) g2 »- 1/2 7, 



-2n-5/6 



2H/6 r 2 (1/2 -n) i 



exp 



2y/6> 



7 F 



In 



1 



(93) 
(94) 

(95) 



We recall that, due to the approximations performed in the previous derivation, this estimate is valid only for 
/i <C /imax- Nevertheless, the argument in the exponential amplifies the transition from bound states to tunneling 
ones for masses /i ~ /i mal , as intuitively expected. An order of magnitude of the minimal coupling constant 7 F leading 
to stable bound states can thus be estimated by requiring that the lowest massive state does not tunnel, 



jUi < 2^ max exp [ -1 -= 



1 



V6 2v / 6#i 



Using the two values (78) and (|S0D, this implies that 

7 F > Hi exp ( 2 



\/6 yfiHx. 

As a numerical application, for the Higgs and gravity parameters used in Fig. H, we get 7 F > 25. 



(96) 



(97) 
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V. NUMERICAL INVESTIGATION 



Numerically, it is simpler and more convenient to solve the first order differential system pQ-pl|). A Runge-Kutta 
integration method was used, on both sides of the brane. In order to suppress the exponential growth, we integrate 
from the turning point, g — I, where the solution begins to propagate freely, toward the brane. In this way, we get 
only U (a, x) near the brane. The radiative solution for g > £, is simply obtained by integrating from the turning 
point toward infinity, with initial conditions determined by the matching with the exponential decreasing solution 
near the wall. The same method is used on the other side of the brane, but this time, by means of the last free 
parameter, we impose the continuity of one bulk spinor component on the brane (U^ say). Generally, the other bulk 
spinor component will be discontinuous at g = 0, as expected from the analytical study since U'(a,0) is generally 
discontinuous. The mass spectrum is thus obtained by requiring the continuity of Un on the brane. 

The bulk spinor components computed this scheme have been plotted for the first massive modes trapped on 
the brane in Fig. pi for 7 P = 100. The lowest mass is numerically found to be /ii ~ 0.2097 F and was estimated 
analytically, from (j78|), to be fii ~ 0.2107 P leading to a precision of 0.5% for the analytical estimate. The second mass 
is numerically found to be fj,2 ~ 0.29l7 F , which has to be compared to its analytical estimate [1% ~ 0.2957 F . Again, 
the precision of this estimate is of about 1%. As predicted from (H), there are n max = 11 massive bound states the 
lightest masses of which are summed up in table |[ On Fig. |[ we plot the last n , = 11 trapped mode; it has a tiny 
radiative component, as expected for a tunneling mode. 




_2 J , , , , , 1 -1e-07 J , 1 

-1.5 -1 -0.5 0.5 1 1.5 -3 -2 -1 

P P 
FIG. 5. [Left] The right (dashed line) and left (solid line) bulk spinor components Wr and Wl, as functions of the dimensionless 
distance g to the brane, for the heaviest massive state (n = 11) and 7 F = 100. [Right] Zoom near the turning point and transition 
to the radiative behaviour which takes place at a finite distance to the brane due to the tunneling of this mode from the brane 
to the bulk. 



n 


1 


2 


3 


4 


5 


6 


11 


[in (numerical) 


0.209 


0.291 


0.353 


0.402 


0.444 


0.480 


0.593 


[i (estimates) 


0.210 


0.295 


0.359 


0.412 


0.458 


0.499 


0.657 


precision (%) 


0.5 


1.3 


1.7 


2.4 


3 


3.8 


9.7 



TABLE I. Comparison of the numerical values and analytical estimates of the first six bound states reduced mass /i, together 
with the heaviest mode, computed for 7 F = 100. 
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In conclusion, the numerics confirm that the approximations of the previous section and our estimates are accurate 
up to 1%-10% (see table @). 

2 i 1 2 , 1 




-0.5 



0.5 



-0.5 



0.5 




-0.5 



0.5 



-0.5 



0.5 




-0.5 



0.5 



-0.5 



0.5 



FIG. 6. The right (dashed curve) and left bulk spinor components, Ub. and Ul, as functions of the dimensionless distance q 
from the brane for the lightest massive bound states. They have been computed for 7 F = 100 and in the Higgs and warp factor 
profiles obtained with parameter of Fig. ^j. The numerical values of the corresponding reduced masses are reported on table [l| 
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VI. DISCUSSION AND CONCLUSIONS 



We shall now discuss the cosmological constraints existing on the kind of model we have been considering here. 
Most of these constraints come from brane models in which the wall structure is replaced by an infinitely thin four 
dimensional layer. As discussed in the previous sections, such an approximation is equivalent, within our framework, 
to asking that the combination a(3 be much larger than unity. In this limit, equivalent to the large (3 limit since, from 
Eq. (|l5|), a/3 ~ we can replace the stress-energy tensor (||) by the effective four dimensional surface distribution 



whose isotropic tension is obtained by integration in the transverse direction to yield 



Too = V\Mv 2 J dgt 



-6<r(g) 
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where the function can be expressed as 



6S 2 -1 



2V 



2/3(H 2 - if 



(98) 



(99) 



(100) 




FIG. 7. Effective four dimensional brane tension, in units of y^\A\r/ 2 , of the domain wall as a function of j3 (dashed line). It 
is clear that in the thin brane limit, £ ~ (5 (thin line). 



Fig. fj] shows the variation of £ as a function of (3. In the limit 3> 1, it is clear that £ ~ (3, so that the brane 
tension behaves as 



Xi] 4 



(101) 



which will be used to derive the relevant cosmological constraints. Note also that the discrepancy between the above 
formula and the actual value of becomes important (more than 100 % error say) for (3 % 0.1, which is already 
rather far from the thin brane limit usually considered. 



A. Investigation of the parameter space 



The model described in this article depends of five parameters, four describing the spacetime and scalar field 
dynamics (G 5 ,A, 77, A) and one concerning the fcrmions (<7 F ). With the domain wall structure assumed, only four of 
these parameters are independent [see Eq. ([l6])]. It is convenient to replace this set of parameters by the three mass 
scales 



1G 



-A 



(i 



in _ 



T l/4 



(102) 



and the dimensionless parameter 7 F . These parameters are subject to a number of constraints, namely 



1. the four dimensional gravitational constant must agree with its observed value G 4 = to~ 2 with m 4 ~ 10 19 GeV. 
Using the expression of the four dimensional Planck mass in terms of the five dimensional analog and of the 
brane tension gives fiE2,E3| 



to 3 ~ to. to 2 



(103) 



2. The brane cosmological constant p 22 2j| 



2A 4 = A + 6n G^T^ , 



(104) 



must also agree with the standard observational bound toa 4 < 10 to 4 . This implies 



2 4 
7T TO^ 



™4 ~ m A m|. 



(105) 



Note that in the limit /3 3> 1, this relation is equivalent to Eq. (101). This means that this condition is readily 



satisfied in the thin brane limit. At this point, it is worth emphasizing that this is precisely the limit in which 
the analytic approximation for fermion masses are the most accurate. 

3. There must not be any deviation of the law of gravity on the brane with respect to the inverse square Newton 
law above 1 millimeter pi|. This implies pa] 



7tl . > 10" 3 eV. 



(106) 



4. Finally, we require the fermion stress-energy tensor to be negligible with the brane stress-energy, so that we 
impose that the mass of the heaviest fermion is smaller than the brane mass scale. By means of Eq. ([so]) , this 
condition reads 



where Hi ends up being function of (3 only by means of Eqs (|T^) and (| 

6/3 



H 



2/3- 



V16/3 + 9 

In the limit (3 ^> 1, Eqs. (|107|) and (|108|) combine to give the constraint on the coupling constant 



7 F < — e 

TOA 



TOA 



(107) 



(108) 



(109) 



It follows from the relations (103) and ( |105| ) that m 2 ^ ~ toaTO4 so that the three mass scales must satisfy 

TO A >10" 3 eV, m^^lTeV, m 5 > 10 6 TeV, (110) 



which, together with Eq. (10S) and ( p5| ) yields 

1 



2y/6 



< 7 F < 10 15 



ITeV 



(111) 



where, as discussed below Eq. (40j), the lower bound is very conservative. 

Now, let us examine the stability of the fermion confinement on the brane and the restriction on their life-time 
imposed by the previous conditions. We require, at least, one massive bound state to have a life-time longer than the 
age of the Universe, i.e., 
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Tl > T u 



10 17 s 



(112) 



Using Eq. (|95|), one roughly gets 



The condition ( |112| ) together with the former constraint ( JllOf ) can be written as 

T(/3, 7f )> 10 30 . 



(113) 



(114) 




FIG. 8. Contour plot of the function T(/3, 7 F ) for T = 10 30 [solid lines] and 10 60 [dashed line], i.e., respectively a particle 
life-time of the order of the age of the Universe, and of the proton life-time lower limit. Note that both top curves are 
indistinguishable due to the exponential behavior of T in 7 F . They turn out to be equivalent to the analytical requirement 
given by Eq. (^). Points above the highest curve and below the lowest curves satisfy the contraint. We also superimpose the 
conservative constraint (^) necessary for the existence of a massive bound state, the allowed region being above the dotted 
line. 



On Fig. |[ we present the contour plot of the dimensionless function T(/3, 7f ), for T = 10 30 and 10 60 , which 
correspond respectively to a particle life-time of the order of the age of the Universe, and of the proton life-time lower 
limit. For /3 > 1, there are in principle two allowed regions, corresponding to strong and weak coupling limits, i.e. 
7f S> 1 and 7f <C 1. However, the lower bound on 7f , which comes from the requirement that fermions are actually 
trapped on the brane, pushes the weak coupling allowed region to very high values of f3, in practice (3 > 7.5 x 10 57 for 
T = 10 30 . Note also that this already rather extreme value is based on the conservative estimate given by Eq. (|35|). 

For P < 1, the weak coupling region completely disappears, while the strong coupling allowed region shrinks 
rapidely: for (3 < 8.5 x 10 -3 , the life-time cannot exceed the age of the Universe because 7f > 10 15 . Considering 
/3 ^ 1 therefore turns out to be the relevant limit if one wishes to have fermionic bound states living on the brane. 



B. Conclusion 



In this article, we have considered fermions coupled to a Higgs field with a domain wall structure in a five dimensional 
anti-de Sitter spacetime. This domain wall can be thought of as a realization of a brane universe. 

After, studying the domain wall configuration, we solved the Dirac equation and showed that there exists massive 
fermionic bound states trapped on the wall. We develop both analytic approximation to compute the mass spectrum 
and the tunneling time. This was compared to a full numerical integration of the dynamical equations that revealed 
the accuracy of our approximation scheme. 

We recover the fact that massive fermions tunnel to the bulk jKj. Investigation in the parameter space shows that, 
for models satisfying the cosmological constraints, the relevant confinement life-time can be much greater that either 
the age of the Universe or the proton life-time. This was made possible by the derivation of the analytic estimate. 
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One of our central result is the derivation of an analytic mass spectrum for fermions trapped on a brane-like four 
dimensional spacetime. In particular, as could have been anticipated fl3|| , it was shown that the allowed masses are 
quantized, with a spectrum varrying, in the strong coupling limit, as ^fn. Such a spectrum is indeed in contradiction 
with experimental measurements of particle masses |26| , which is not surprising given the simplicity of the model. It 
however opens the possibility to build more realistic theories in which mass quantization would stem naturally from 
extra dimensions. 
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